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The methods of non-equilibrium thermodynamics of systems with an interface have been applied 
to the study of transport processes in semiconductor junctions. A complete phenomenological model 
for drift-diffusion processes in a junction has been derived, which includes, from first principles, both 
surface equations and boundary conditions, together with the usual drift-diffusion equations for the 
bulks. In this way a self-consistent characterisation of the whole system, bulks and interface, has 
been obtained in a common framework. The completeness of the model has been shown and a simple 
application to metal-semiconductor junctions developed. 



I. INTRODUCTION 

The phenomenological approach to semiconductor de- 
vice modelling has proven to be very powerful and 
successful during the last decades. By means of a 
relatively simple set of equations (i.e., balance equa- 
tions, constitutive relations, Maxwell equations), numer- 
ous devices have been modelled jl). In spite of this suc- 
cess some aspects of the phenomenological models still 
remain not so well developed, namely the boundary con- 
ditions necessary to solve the set of phenomenological 
equations and the appropriate equations to describe the 
surface processes phenomenologically. These are, by no 
means, secondary problems in semiconductor devices, 
since the behaviour of most of the devices are critically 
affected by them. For instance, any device including a 
junction needs the appropriate boundary conditions, or 
discontinuity relations, to be imposed at the junction, as 
well as, the proper modelling of the superficial phenom- 
ena. In the term junction we also include the contacts 
used on the devices, whose effects may sometimes be rel- 
evant, for instance, in the development of semiconductor 
instabilities 

Usually, these problems are discussed more or less in- 
tuitively, by invoking sometimes unclear physical argu- 
ments. We may say that, at this moment, in the current 
literature, a general phenomenological theory for bound- 
ary conditions and surface effects is lacking, although 
some valuable work has been published in this area ||. 
It seems, therefore, desirable to propose phenomenolog- 
ical models that incorporate both boundary conditions 
and surface equations in a fundamental and general way. 
With such models, a self-consistent characterisation of in- 
terfaces in semiconductor systems would then be able to 
be reached. Therefore, for a given junction, the models 
would provide the corresponding surface equations and 
boundary conditions appropriate for it simply by adjust- 



ing a few surface parameters. Clearly, from this point of 
view, semiconductor device modelling turns out to be a 
much easier task. Moreover, only with this kind of ap- 
proach, the presently accepted models will be ultimately 
justified and, if necessary, improved in a systematic and 
general way. 

As shown in ||, these kinds of models may be con- 
structed by means of non-equilibrium thermodynamic ar- 
guments. In this paper we will also make use of such an 
approach, but following a rather different formulation; 
namely, the one proposed in J(J, in the context of fluid 
systems (see review and references quoted therein). 
Non-equilibrium thermodynamic derivation of transport 
models have been widely used in the past for bulk sys- 
tems (J], usually in the context of fluid systems, but also 
in the context of semiconductor systems || . For systems 
with interfaces, much less work has been done, although 
a rather complete theoretical treatment of this problem 
has been already published JjJ . 

In this paper we will present an application of such 
a formalism, to the derivation of a transport model for 
drift-diffusion processes in a junction. As required above, 
this model will incorporate, apart from the usual drift- 
diffusion equations for the bulk systems forming the junc- 
tion, the appropriate boundary conditions and surface 
equations to respectively describe the coupling between 
the bulk systems and the surface processes. As men- 
tioned, the entire model will be obtained in a systematic 
way, starting from very general considerations about the 
non-equilibrium thermodynamics behaviour of the sys- 
tem. For this reason, its generalisation to include other 
irreversible processes may be obtained following the same 
procedure as the one presented here. 

This paper is distributed as follows: in Section II 
we will discuss in detail the macroscopic description of 
the junction and the semiconductor model we will use 
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throughout the paper. In Section III we will present and 
derive the phenomenological model, and consider some 
particularly interesting cases. The completeness of the 
model will be proven in Section IV, together with a simple 
application to metal-semiconductor junctions. Finally, in 
Section V we will summarize our main results. 



II. MACROSCOPIC DESCRIPTION OF A 
JUNCTION 



From a macroscopic point of view, a junction may be 
defined as the very narrow region that exists between 
two semiconductors in contact, and which exhibits pecu- 
liar macroscopic behaviour reflected in the fact that the 
equations of state, transport coefficients, etc. are differ- 
ent from those of the bulks. It is clear that due to the 
possibility of jumps in the relevant quantities, a descrip- 
tion of the system in terms of variables varying contin- 
uously from one region to the other is not appropriate. 
Therefore, it is necessary to find the jump relationships 
coupling the different parts, in terms of the interfacial 
properties. 

To this end we will introduce in our system a geo- 
metrical surface inside the interfacial region that divides 
the system in two separate parts, referred to as the bulk 
parts, which will be denoted by — and +. We assume 
that each of these bulk parts behaves macroscopically as 
in the absence of interface. Of course this behaviour will 
be different from the real one inside the interfacial region. 
To take this fact into account one introduces surface vari- 
ables, denoted by the superscript s, defined only at the 
geometrical surface. These variables are the appropriate 
to describe the interfacial region. 

The precise definition of these surface variables runs 
as follows. Let 9 + (6~) be Heaviside's functions whose 
values are 1 (0) for points located on the right (left) of 
the surface and (1) for the opposite case. Let d be an 
unspecified field and d be the corresponding fields in the 
bulks, referred to as the extrapoled fields. As pointed out 
before, these fields differ from the real ones in the inter- 
facial region. To account for this difference we introduce 
excess fields d ex as follows H 

d ex {r) = d{r) - d+(r)6+(r) - d~ (f)0~ (f) (I) 

These excess fields have, by definition, two important 
properties: 1) They vanish away from the interfacial re- 
gion, because there the effects of the interface die out, 
and consequently the real fields and the bulk ones co- 
incide. Thus the excess fields are only defined in the 
interfacial region. 2) They take small values when inter- 



facial effects are small, that is, when the real fields do not 
depart significantly from the extrapoled ones. 

The surface variables are defined as averages of the 
excess fields. For a planar surface at rest identified with 
the plane x = 0, we can write JT(J 

d s (y,z) = J dxd ex (f) (2) 

where the integration may be extended outside of the in- 
terfacial region, for the excess fields vanish very quickly 
there. These fields constitute the relevant quantities in 
the description of continuum systems in the presence of 
interfaces, as was first pointed out in |@). 

Note that, with this characterisation, the system is 
seen as composed by three sub-systems: two bulks, with 
unchanged properties just until the geometrical interface 
and described by means of the extrapoled variables, and 
an interface, accounting for the properties of the interfa- 
cial region through the surface variables. 

As noted before, in this paper we will follow a ther- 
modynamic approach. This forces us to give an a priori 
thermodynamic characterisation of the systems involved. 
Here we will assume that the bulk sub-systems belong to 
the class of the multicomponent electrically conducting 
polarizable fluids whose thermodynamic variables are: 
the densities of internal energy, u, and entropy, s, the 
number density of each component, Uk, and the electric 
polarisation P (magnetic polarisation, M, may also be in- 
cluded). The surface sub-system will belong to the same 
class, but neglecting any polarisation phenomena. The 
bulk system densities will be defined per unit of volume 
while at the surface per unit of area. 

This semiconductor fluid-like model we will use 
throughout this paper, is very related to the one pro- 
posed in H in which a semiconductor may be considered 
as the superposition of three continuum media, namely 
the crystal lattice, the electrons and the holes. The first 
one is, in general, positively or negatively charged de- 
pending on the sign of the ionised impurities. It will be 
considered at rest with fixed mass and impurity densities, 
not including for example any aging phenomena. This 
sub-system is mainly responsible for the polarisation ef- 
fects. Electrons and holes will be considered charged flu- 
ids. For the sake of simplicity, we will assume that each 
component has the same charge and mass on each side 
of the junction. This assumption considerably simplifies 
our analysis in the sense that we are reducing the number 
of constituents in our description. 

It is worth nothing that although the formalism we 
present applies to semiconductor systems, it can also be 
used to study other kinds of materials, like metals or in- 
sulators, because from the thermodynamic point of view, 
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they have similar characteristics. It is only when the 
state equations are specified that the differences become 
substantial. 



III. NON-EQUILIBRIUM THERMODYNAMICS 
MODEL FOR A JUNCTION 



Once the thermodynamic description of the junction 
has been established, we can now proceed to investigate 
the non-equilibrium behaviour of the whole system by 
applying the usual methods of non-equilibrium thermo- 
dynamics B . Our main purpose in this section will be to 
obtain the constitutive relationships for the irreversible 
processes in the bulk and at the interface together with 
boundary conditions coupling the different parts of the 
system. As noted before, a similar approach was previ- 
ously proposed in |j) in the context of fluid systems. 

Although the calculations may be worked out in the 
very general case, by including all possible irreversible 
processes as in ||, we will consider, for the sake of 
simplicity and clarity, only the drift and diffusion pro- 
cesses, neglecting the recombination or generation pro- 
cesses, magnetic effects, thermal processes or relaxation 
polarisation processes. The last three may be intro- 
duced in a more or less straightforward way by follow- 
ing the procedure presented in this paper. Generation- 
recombination processes deserve special treatment in the 
framework of non-equilibrium thermodynamics, due to 
the intrinsic non-linear nature of the constitutive rela- 
tions 101. 



Poisson's equation, supplemented by the corresponding 
equations of state. They have been found to be 
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dt 
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dD ± _ ± 
dx 9 
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(3) 
(4) 
(5) 
(6) 



where n k is the particle number density with the index 
k referring to the subcontinua electrons (k = 1), holes 
(k = 2) and lattice (& = 3), = J2k=i Qk4 ^ s the total 
density charge, with q k the charge of the fc-th compo- 
nent, is the electric displacement, JjT is the particle 
number density current, is the diffusion coefficient 
of fc-th component, and jl k the corresponding mobility. 
These two last quantities are related through Einstein's 
relation 
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k = 1,2 



(7) 



where e is the (positive) electron charge. These bulk 
equations constitute the usual drift-diffusion model for a 
semiconductor when generation-recombination processes 
are neglected (Eq. (Q) simply states that the aging 
phenomena are neglected) . Eq. (g) can be rewritten in a 
more appealing form by introducing the electrochemical 
potential for the fcth component as 
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(8) 



Moreover, we will also consider a one-dimensional sys- 
tem in which significant variations of the relevant quan- 
tities used in its description occur only in a direction 
perpendicular to the interface. Vectorial fields will then 
have components along the normal to the junction, and 
the vectorial notation will be suppressed. 



with the electric potential defined as usual by E = 
—dV /dx. For the constitutive relations, one then ob- 
tains 



-44' 



dx 



1,2 



(9) 



A. Basic equations of the model 



This subsection is devoted to the introduction of the 
complete set of equations that constitutes the model, 
leaving its rigorous derivation for next subsection. As 
mentioned in the introduction, the equations of the model 
may be divided in three categories making reference to 
bulk, surface, and boundary conditions, respectively: 

i) Bulk equations. The equations describing the bulks 
consist of balance equations, constitutive relations and 



which shows that, for drift-diffusion processes, the elec- 
trochemical potentials (usually called quasi-Fermi levels) 
play an important role, as we will see later when bound- 
ary conditions are expressed in a natural way in terms of 
these potentials. 

The appropriate equation of state for the bulk systems 
corresponding to the model presented are 



4 
D ± 



4(4) 

e ± E ± 



k = 1,2 



(10) 
(11) 



where are the bulk dielectric permitivities, assumed 
to be independent of n k . For the case of non-degenerate 
semiconductors one has for the first state equations 
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Mfc = M fc ,o + k B T In 



k = 1,2 



(12) 



where /i^ is a chemical potential of reference, with 
being the effective density of states in the conduction 
(k = 1) and valence (fc = 2) band. In that case Ein- 
stein's relation reduces to its usual expression [Q , 
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k = 1,2. 



(13) 



1/2 (A + + A + ), and where all the bulk magnitudes are 
evaluated at the interface, that is, at x — 0. 

These equations can be formulated in a more com- 
pact way, by introducing a simplified notation for the 
differences of electrochemical potentials, X k = (f> k — <^|, 
X k = <p A k — (p k , called, from now on, surface thermo- 
dynamic forces, and for the fields corresponding to the 
same quantity but for different constituents, for instance, 
J ± = (Ji, J^) where 1 and 2 refer to electrons and holes 
respectively. The scalar product of two such vectors will 
then read J + • X + = J k X k . In this notation, 

Eq. (|l^ ) can be written simply as 



ii) Surface equations. As for the bulk, surface equations 
would consist in general of balance equations, constitu- 
tive relations and state equations. For the case of a one- 
dimensional system no surface constitutive relations are 
needed, because no transport process takes place along 
the surface. In that case, we then simply have, 



-dT + [Jk] - 
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dt 







(14) 
(15) 



where n s k is the surface density number of the kth com- 
ponent, and where we have defined [A]- = A + — A~, 
with A^ being evaluated at the interface. Again Eq. Jl5| ) 
states that no aging phenomena are allowed, now at the 
interface. Moreover, the surface state equations are in 
that case simply 



/4 = 



1,2. 



(16) 



Remember that we are assuming that the interface is 
non-polarizable and hence no surface state equation for 
the surface polarisation is necessary (see appendix) 

iii) Boundary conditions. For the boundary conditions, 
we have to distinguish between the ones corresponding to 
Poisson's equation and the ones that correspond to the 
drift-diffusion processes. The former ones may be written 
as 



[D]- = q S 
[V]-=0 



(17) 
(18) 



where q s = X)fc=i Qkn s k is the net surface charge, whereas 
for the latter we have obtained 



J. 



where 



F t(i^j ~ ^ ~ 4>J}) ; k,j = 1,2 (19) 



M| + is the electrochemical poten- 



tial of the fc-th component at the surface, with [A]- 



F ± (X+,X- 



(20) 



The boundary conditions for the drift-diffusion processes, 
Eq. ( |T9| ) or Eq. (p0|), deserve some comments. First of all, 
it should be noted that they are non-equilibrium bound- 
ary conditions, that is, they relate the discontinuities 
in the relevant quantities in the drift-diffusion process, 
namely the electrochemical potentials to the currents. 
Moreover, although they have been formulated in terms 
of currents and forces (see above), they may also depend 
on the surface thermodynamic variables (see next sub- 
section). Finally, the functions are subject to certain 
thermodynamic restrictions, namely 



^(0,0) = 

X + ■ F+ + X- ■ F~ > 



dFS 



dX 



dFf 



dxt 



eq 



fc,i = i,2 

a, (3 = +,- 



(21) 
(22) 

(23) 



coming from the equilibrium condition, the second prin- 
ciple of thermodynamics, and the Onsager's symmetries, 
respectively. Note that we have n(n — l)/2 symmetry 
relations, where n is the number surface thermodynamic 
forces (in the general case n= 4). 

It should be noted that for the model presented above, 
the total electric density current is a continuous field 
through the interface. Indeed, from Eqs. ( |T4| ) and (|T5|), 
the total surface charge equation for a one-dimensional 
system is found to be 



J^i + [J e ]_ = 



(24) 



where J^r = Y], quJ^ , is the total electric conducting 
density current. Upon substituting Eq. ( p~7[ ) in Eq. (|2J 
we obtain 
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[7]_ = 



where 



(25) 



(26) 



is the total electric density current which can be shown 
from the semiconductor bulk equations, to only be a func- 
tion of time (Ampere's law). This fact, and Eq. (25) 
makes the use of a superindex for this quantity irrele- 
vant. 



where Lif^ is the matrix of phenomenological coefficients 
which according to Onsager's principle is symmetric and 
according to second law positive definite [I] , and only de- 
pends on the bulk thermodynamic variables, which in our 
case are, for instance, T^, rvt, £ ,± . To arrive at the form 
given in Eq. (||) we will suppose that the cross effects 
are negligibly small, Lf 2 = L^i ~ and that the phe- 
nomenological coefficients and the chemical potentials do 
not depend on E ± . We will also assume that the chemi- 
cal potentials depend only on the corresponding density, 
that is //£ = [if (nf ). Then defining the diffusion coeffi- 
cient and the mobility for the electrons and holes as 



B. Entropy production and constitutive equations 



this subsection is to justify the equa- 



Our purpose in this subsection is to ju 
ans proposed in the previous subsection 



tions pro 

i) For the bulk equations we only have to derive the 
constitutive relations, Eq. (|J), and Einstein's relation 
Eq. (0), because Eqs. (0) and (Q) are simply continuity 
equations and Eqs. (|10|), (fill), appropriate state equa- 
tions. More general equations of state would be possible 
a priori, but we will conclude from our analysis that these 
are the consistent ones with the remaining equations of 
the model. To derive the constitutive relations, we will 
analyse the bulk entropy productions using the methods 
of non-equilibrium thermodynamics Q] . In the case con- 
cerning us, that is when only drift-diffusion processes are 
considered, the entropy production reduces to (see ap- 
pendix) 



f 



T ^ 

k=l 



d/M k 



Ox 



(27) 



where the quantities Jjjr are diffusion currents defined in 
the appendix. After applying Prigogine's theorem j|, for 
which we can redefine the diffusive fluxes with respect to 
any desired velocity (in our case we will take the reference 
frame fixed to the lattice continuum which is assumed to 
be at rest), Eq. (p7|) transforms into 



T J k 

k=l 



dx 



(28) 



From it we can infer the phenomenological constitutive 
relations. For bulk processes they are well represented 
through linear relations, that is, 



f 



T ^— ' 



fcj V dx 



k = l,2 (29) 
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k — 



T Or 



L 



kk ' : 

T n 



k = 1,2 



k = 1,2. 



(30) 
(31) 



Eq. (g9J) transforms into the drift-diffusion relation given 
in Eq. (§. 

Einstein's relation, Eq. (Q), follows directly from the 
previous definitions of the diffusion and mobility coeffi- 
cients. On the other hand, taking into account the def- 
inition of the electrochemical potential given in Eq. (||), 
we may rewrite Eq. (|29|) as 



1 2 



dx 



k = 1,2 



(32) 



which after neglecting cross effects and considering the 
definition of the mobility given through Eq. (|3l]), trans- 
forms into Eq. (m. 



ii) The surface equations, Eqs. ( |i4| ) and ( |l5| ) are the 
corresponding one-dimensional surface continuity equa- 
tions and they have been derived in the appendix . As 
before, 



Eq. (10) corresponds to particular forms of the 
surface state equations appropriate to the model. 



iii) Finally, let us focus on the boundary conditions. 
Eqs. (|l7| ) and (fl8|) are the usual discontinuity relation 
for Poisson's equation and they do not deserve further 
comments (see the appendix for the case in which sur- 
face polarisation is included in the formulation and (Tc| l 
for a general treatment of Maxwell's equations). 

To derive the remaining boundary conditions, Eq. (|Ts|), 
we proceed, as we did in the bulk, by analysing the sur- 
face entropy production. It is shown in the appendix that 
for the case of drift-diffusion processes, it may be written 
as 
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°i=-^iy* fa -&)]-> (33) 
fc=i 

Using the definition of the bracket [•]_ we may rewrite 
Eq. ( p3| ) as 

2 2 

^ = ~E Jk(t4 f4) \ E - Mfe ) • (34) 

fc=i fe=i 

Finally, by using Eq. (|l^) and the definition of the surface 
and bulk electrochemical potentials given in the previous 
subsection, we arrive at our final expression for the sur- 
face entropy production, 

2 2 

< = ~ E -tut - - f E -TO* - ^) ( 35 ) 
fc=i fc=i 

= -Ij+.X+-ij--X-. (36) 



C. Particular cases 

1. Steady conditions 

Under steady conditions, the model presented in sub- 
section III. A, and in particular the boundary condi- 
tions, simplify considerably, and they can be written in 
a more useful form for practical applications. Indeed, 
from Eq. @J) we conclude that no discontinuity in the 
conducting currents occurs across the interface, that is , 
,/ + = J~ at the interface. Hence, after identifying the 
terms of the right hand sides in Eq. ( |20| ) and solving for 
one of the two forces, we can write the boundary condi- 
tions in the form 

X- = G+(A+) (37) 
J+ = -IiT+(X+) (38) 



This last expression enables us to identify the surface cur- 
rents, , and their conjugated surface thermodynamic 
forces, X in a clear way; hence, we can propose the phe- 
nomenological relations between them. As we will see 
in the application, for interfaces, general kinds of phe- 
nomenological relations need to be used, and for that 
reason we have written them like in Eq. ( |20| ) , where the 
dependence of the functions F , on the surface thermo- 
dynamic variables, T s , n s kl is implicit. 

Equilibrium and non-equilibrium thermodynamic con- 
siderations impose some restrictions on the functions F^ , 
namely: 

a) At equilibrium, all thermodynamic fluxes and 
forces must vanish, which leads to the requirement 
that F^ have to vanish in equilibrium, that is, 
F ± (0,0) = 0. 

b) They have to be defined in such a way as to make 
surface entropy production positively definite, in 
accordance with the second law of thermodynam- 
ics, that is, X+ ■ F+ + X- ■ F- >0 (see Eqs. @ 
and @ ). 

c) Near equilibrium, in which linear phenomenologi- 
cal relations hold, the matrix of phenomenological 
coefficients satisfy Onsager's reciprocal relations. 
This matrix is precisely the Jacobian matrix for the 
functions F^ evaluated at equilibrium and, conse- 
quently, it must have the corresponding symmetry 
properties. 

These three requirements are, in fact, the conditions ex- 
pressed through Eqs. @- @. 



where G + is a function satisfying the relation 

F + (X+,G + (X+)) =F-(X + ,G+(X + )) (39) 



and H + is given by 

H+{X+) = F + {X+, G+{X+)) . 



(40) 



It can be seen that the general restrictions mentioned 
previously require that G + (0) = # + (0) = and that 



{X+ + G+{X+)) ■ H+(X+)>0 



(41) 



for all values of X + . 



This last relationship is more naturally expressed by 
using the discontinuities in the bulk electrochemical po- 
tentials, [</>]_, as the independent variables. In this case, 

by taking into account that [(f)] _ = X + + X~, we may 
rewrite Eqs ( |37| ) and (|3^) as 



J+ -- 
X+ 



1 

F 2 ( 



where F 2 satisfies 

G + (F 2 M_))+F 2 (14>}_) = {4,}_ 

and F\ 



(42) 
(43) 



(44) 



G 



(45) 



IV. APPLICATION OF THE MODEL 



Now the thermodynamic restrictions are simply 

F 2 {0) =F 1 (Q) =0 (46) 
m_-F 1 M_)>0. (47) 

The steady state boundary conditions expressed through 
Eqs. (f|2|) and ([43|), have the appropriate form to be used 
in practical calculations, as we will show in future works. 



2. No surface states 



In this section we will first show, in a more or less 
heuristic way, that the model presented along this paper 
is complete (see below). For simplicity's sake we will fo- 
cus on unipolar systems, but similar arguments can be 
applied to bipolar systems. As an example of application, 
we will show that the thermionic emission-diffusion the- 
ory for metal-semiconductor (M-S) abrupt junctions [EJ, 
is a thermodynamically consistent model, because it can 
be fitted into the framework of the present model. Simi- 
lar arguments may be used to analyze more complicated 
situations, in which not only thermodynamic consistency 
is expected to be proven, but also a generalization of the 
corresponding models will be derived. 



Our purpose in this subsection is to present the sim- 
plified version of the model corresponding to the case 
in which the surface number densities vanish, that is , 
when n\ = for i = 1,2,3. This case may be related 
to the case in which there are not surface, or interfa- 
cial, states. The first simplification comes from the fact 
that now we obviously, do not need a surface state equa- 
tion like Eq. (|l^). Moreover, from the surface continuity 
equations, Eq. (jlj), we can see that the number den- 
sity currents are continuous across the interface, that is, 
[</]_ = 0. Finally, it can be shown that the entropy pro- 
duction corresponding to this situation can be obtained 
from Eq. ( fj3| ) by simply imposing [J]_ =0 on it. One 
then arrives at, 



< = —[$]- -J' 



(48) 



From this expression we see that, in this case, we obtain 
a single discontinuity relationship, that can be written as 



(49) 



J + = ~M{[$}.) 



instead of the two obtained in the general case Eq (J2C 
or Eqs (|42] ) and (fl|) in the steady case. As before, M is 
subject to some thermodynamic restrictions, namely 



M (0) = 
$_-M($_)>G 



dMi 



d[<h 



dM 2 



d[4>i 



(50) 
(51) 

(52) 



coming from the equilibrium condition, the second prin- 
ciple of thermodynamics and the Onsager's symmetries. 



A. Completeness of the model 



What will be proven in this subsection is that if the un- 
specified phenomenological equations, Eqs. (|l6|), and (|2(i| ) 
in the general case or Eq. (|49[) in the non-surface state 
case, are given explicitly and the physical conditions at 
infinity are specified, then we can find the solution of 
the model. By physical conditions at infinity, we mean 
those conditions that one assumes the solution must sat- 
isfy far away from the interface. For instance, for infinite 
systems, one assumes that the interfacial effects die out 
at infinity, and hence that the system is locally neutral 



there, that is, — > n$ and -q£- — > for x — > ±oo, 
where we have assumed the single carriers to be elec- 
trons (subindex e) and the donor density (subindex 0) to 
be uniform. For finite systems, one assumes that, at a 
certain distance away from the interface, the interfacial 
effects have nearly died out, and, therefore, one searches 
for an approximate solution that satisfies this condition. 

In any situation, our model provides the correspond- 
ing solutions. To see this, let x^ be, for the moment, two 
unspecified points, at which the neutrality conditions are 
satisfied within a certain degree of approximation, that 



n e (x , t) 
dnf(x^t) 
dx 



At 



(53) 
(54) 



where A^ 2 take small values ( for infinite systems these 
precise values are not important, for at the end of the 
argument one tends them to zero). Note that by taking 
Eqs. (53) and ( |54| ) as boundary conditions, one can solve 
the semiconductor bulk equations, Eqs. (|3|), (||), ([|). A 
simple way to do this is to note that that these equations 
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can be transformed into a single equation for the electric 
field, see for instance El] ], 



~dt 



-Dz 



d 2 E ± 

dx 2 



P ± E ± 



(dE A 



\ dx 



I 

7± 



(55) 



where / is the total electric density current, and that 
Eqs. (HH) and (|54| ) can be rewritten as, 

0£±(4,t) 



dx 



£±(4,*) = — jf{xt,t) 

/T->- 



(56) 
(57) 



where Eqs 
and ©, c 



(I) 

.± _ 



and (c 

e/i ± n 



) have been used. In Eqs. 
are the bulk conductivities, 



|5£ 

and 



A^ 2 , two very small quantities, related to A^. 



Once one has solved for the semiconductor bulk equa- 
tions, (usually numerically), one can compute some use- 
ful quantities like: Jt(0,t), nf(0,t), /J,f(0,t), V^O,*) - 
v± i x oi t )i 4(°>*) - 4>t{ x Q^)^ ■■■ b y simply using the 
bulk equations again. Note that all these quantities will 
obviously be dependent on xj, I(t) and Af 2 . If surface 
states are present, one can then compute the correspond- 
ing expressions for the surface quantities, q s (t) and jJ. & e {t) 



as functions of 



/ and Af 2 . 



follows: from Poisson's equation, Eq 
continuity equation, Eq. 



cqua 



This can be done as 
(||), and the dis- 

one has 



E+ (x+ t) - E~ (x Q , t) = q- (t) + q+ (t) + q s (t) (58) 



where we have defined 



±- 



(nj - nf)dx. 



(59) 



From Eq. ( |58| ) one obtains q s (t), for it is the only un- 
known quantity in this expression. By using the surface 
state equation, Eq. (16), one can then derive the corre- 
sponding expression for n s e (t) (we assume that the surface 
donor density, rip, as well as the bulk ones, , is known). 

All the previous results suffice to derive the expres- 
sions for X^(0,t) and J^(0,t) in terms of xj, / and 
Af 2 . By substituting them into the discontinuity rela- 
tionships, Eq. (|2fj|), one then obtains a system of two 
equations, which allows one to derive the expressions for 
x^ as functions of / and A^ 2 . At this point one is then 
able to eliminate the dependence of all the previous re- 
sults on x^. As a result, all these magnitudes become 
only functions of I and A^ 2 . 

For an infinite system, one then tends A^ to zero, 
(and consequently x will tend to ±oo). The resulting 



expressions, which are now only functions of/, constitute 
the searched solution of the model. 

For a finite system, one assumes that the results found 
with finite values for A^ 2 ) are the (approximate) solu- 
tions for the active region of the device, and that the 
remaining parts of the system behave as in the absence 
of the interface. We then arrive at the desired solution 
in this case. In some cases, for instance for short devices, 
this is not a good approach, and a more elaborated one 
must be carried out JljJ| . In any case, it can be shown 
that the correct implementation of the model gives the 
corresponding solution. 

For the particular case in which there are no surface 
states present, the same derivation holds, by simply not- 
ing that in this case one has q 3 (t) = 0. Therefore, 
Eq. (|5^) constitutes the first equation to determine x^, 
and the second one is the single discontinuity relation for 
that case, Eq. (49). 



Thermionic emission-diffusion theory for M-S 
junctions 



The thermionic emission-diffusion theory for abrupt 
M-S junctions without interface states is a currently ac- 
cepted model to describe the rectifying properties of M- 
S junctions when the main limiting current mechanisms 
are the thermionic emission of carriers over the junction 
barrier and the diffusion of the carriers inside the semi- 
conductor [Q, [^2|. Our purpose in this subsection is to 
show that this model can fit into the framework of the 
model presented throughout this paper, and that, as a 
consequence, it is a thermodynamically consistent model. 

As in many cases, this model assumes drift-diffusion 
transport equations for the bulks Q], which are equiva- 
lent to the ones presented in section III. On the other 
hand, it also assumes the following boundary condition 
to describe the thermionic emission processes, |13 (see 
also §) 



eJe = V R (n m - n+) 



where Vr is a positive constant, and n T 



(60) 



This boundary condition is a thermodynamically consis- 
tent, in the sense of the present paper. Indeed, by using 
Eq. (Jl^) it can be rewritten in the form of a discontinuity 
relation, 



e Jg = Ire 



d0e 



- 1 



(61) 



where It — VrNq. Furthermore, this discontinuity rela- 
tion clearly satisfies the thermodynamic restrictions de- 
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rived for that case, Eqs. ( |50| ) and ([[l]) (for unipolar sys- 
tems without interface states there are no symmetry re- 
lationships). It is worth noting that, although Eq. (5C) 
is expected to be satisfied for all the models, Eq. (51) 
might not, for its origin is purely thermodynamic, and 
it has not been included, as far as we know, in any phe- 
nomenological model. Moreover, it is should be pointed 
out that the constitutive relation given through Eq. j6l| ) 
is effectively non-linear (as we have assumed to be), in 
the sense that it does not relate the current and the force 
linearly. However, this does not constitute a singular case 
in non-equilibrium thermodynamics for there are differ- 
ent irreversible processes (chemical reactions Q , adsorp- 
tion processes |^l| ,...) which display such non- linear be- 
haviours. 



Similar analysis may be carried out for more com- 
plicated models, for instance those including interface 
states. By fitting those models into the framework of our 
model their thermodynamic consistency may be estab- 
lished and, what is more important, their generalization 
eventually derived. These facts illustrate the power of 
the application of the thermodynamic methods to semi- 
conductor systems. 



V. CONCLUSIONS 



In this paper we have presented a fully macroscopic, 
and thermodynamically consistent, derivation of a phe- 
nomenological one-dimensional model for drift-diffusion 
processes in a junction. The model, which consists of 
the usual transport equations for the bulk systems that 
conform the junction, also incorporates the appropriate 
boundary conditions or discontinuity relationships, and 
the equations describing the surface phenomena. These 
last two ingredients are essential from a fundamental 
point of view, and they allow a complete modelling of 
semiconductor devices, in which junctions play a so rele- 
vant role. 

The derivation of the model has been completely based 
on thermodynamic arguments. A first step has con- 
sisted of thermodynamically characterising the junction, 
by means of the introduction of the surface and (ex- 
trapolated) bulk variables. Therefore, the methods of 
non-equilibrium thermodynamics for systems with in- 
terfaces |?J, have been applied giving rise directly to 
the above-mentioned model. It is worth noting that 
this derivation clearly ensures the thermodynamic con- 
sistency of the whole model. Of special interest are the 
thermodynamic restrictions affecting the drift-diffusion 
boundary conditions, Eqs. (pl|)- (p3|), for they must be 
satisfied for any proposed boundary condition. 



As in any phenomenological model, explicit expres- 
sions for some functions or parameters are not given by 
the theory itself. Only when the model is applied to 
several different systems, or when the appropriate mi- 
croscopic models are proposed, information about these 
unspecified parts may be obtained. It should be pointed 
out, that, in any case, these results must satisfy the gen- 
eral restrictions imposed by the thermodynamic deriva- 
tion. 

We have shown, in the particular case of unipolar sys- 
tems, that the model is complete, in the sense that, once 
the physical conditions at infinity are specified, it com- 
pletely characterizes the system. A similar derivation can 
be developed for bipolar systems. 

As a simple application, we have analysed the 
thermionic emission-diffusion theory for one-dimensional 
abrupt M-S junctions without interface states and have 
shown that this theory is thermodynamically consistent, 
in the sense of the present paper, by showing that it fits 
well into the framework of the present model. In ad- 
dition, we have also justified with this example the use 
of non-linear constitutive relations in the formulation of 
the discontinuity relations. Similar procedures can be 
applied to more complicated situations, allowing one to 
extract very useful information on the semiconductor sys- 
tems. 

Finally, it should be noted that the drift-diffusion 
model presented here may be generalised in a more 
o less systematic way, by including other irreversible 
processes or by considering three-dimensional systems, 
which would allow us to study, for instance, inhomo- 
geneous junctions ||l5|| , fusfl . Among these generalisa- 
tions, it is especially important to include recombination- 
generation processes or tunnelling processes in the model 
because their relevance to semiconductor devices is well 
known. The treatment of recombination-generation pro- 
cesses in the framework of non-equilibrium thermody- 
namics will be considered in future papers. On the other 
hand, the non-local nature of the tunnelling processes 
makes their implementation in this (local) framework 
rather complicated, although some local macroscopic for- 
mulations are present in the current literature pl| , p9| . 
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[Jdl- 



(A4) 



APPENDIX A: ON NON-EQUILIBRIUM 
THERMODYNAMICS OF SYSTEMS WITH 
INTERFACES 



Our purpose in this appendix is to establish the non- 
equilibrium thermodynamics basis of multicomponcnt 
electrically conducting polarizable systems divided by an 
interface by following the procedure indicated in (9) . Cru- 
cial in this analysis is the derivation of the surface and 
bulk entropy productions accounting for the irreversible 
processes occurring in the system. The phenomenolog- 
ical equations derived from these quantities correspond 
to the constitutive equations for the interfacial and bulk 
currents and to the boundary conditions. 

First of all, we will indicate the way to derive balance 
equations at interfaces 0. To this purpose, let us focus 
on an unspecified quantity d defined in the whole system. 
Its evolution is dictated by the balance equation 



d 



divJ d = a d 



(Al) 



where Jd and ad are the corresponding current and pro- 
duction, and where the subindex V means that this den- 
sity is per unit of volume. Now, using the decomposition 
of the fields in terms of the bulk fields and excess fields 
given through Eq. (Ry we arrive at 



c) 



—d+ + divJj - a+ 

d 7 j • 7- 

Ql d v + dwJ d -°d 



— d ex + divJ d ,ex + n ■ [Jd] - 5 s - a d ,, 
at 



(A2) 



where the subindex A refers to the fact that the density is 
per unit of area, and where only the parallel component 
of the current appears because the integration of the nor- 
mal component vanishes. Usually one defines a surface 
flux with no normal component and calls this restriction 
a transversality condition. Further details may be seen 
in§. 

We will now apply the former scheme in order to ob- 
tain the set of balance equations describing the evolution 
of the relevant quantities in the system. For the balance 
of mass of the k-th component one has 



0_ 

dt 



Pk + divJ f 



1,2,3 



(A5) 



d 



P s k + divJ s u +n-[J k ]-=<j s k 



1,2,3 (A6) 



where J h 



and are the diffusion currents and a k 



and cr| the production rates. From these equations one 
may derive the balance equations for the total mass 
P ± = ELiPfc' P S = ELiPl and total char ge, Qv = 
ELi z kpf, <1a = ELi z kPk> wit h z k being the charge- 
mass ratio of the k-th component, 



—p +divJ p 



dl P 



divj 



p,\\ 



[J P ]- = 



— q$ + divff = 



divJ*u + n ■ [Jj_ =0 



(A7) 



(A8) 



(A9) 



(A10) 



where 5 s is Dirac's delta function coming from the deriva- 
tives of Heaviside's functions 0, n is the unit normal 
vector pointing from — to + and [A]- = A + (0,y, z) — 
A" (0, y, z) (we are assuming that the interface is at rest 
and is identified with the plane x = 0). In the bulk, ef^ 
obeys the balance equations 



divj. 



(A3) 



As a consequence, the two first lines in Eq. ( |A2[ ) vanish. 
By averaging the remaining third term as well as using 
the definition of the surface fields given in Eq. (0) we 
arrive at the surface balance equation 



where we have taken into account the reaction con- 
servation relations Efc=i a k = ^> Efc=i z k <jZ k = 0' 
Efc=i °fc — 0, Efc=i z k (J k — 0, ancl we have defined the 
total mass currents as — Efe=i Jp \\ — Efe=i ^k \\ 
and the charge conducting currents as = Efc=i z kJk> 

= Efe=l Z kJk,\\- 

In the same way for the balance of momentum one has 
d 



dt 



9$ 



divP* 



fv 



0_ 

at 



? A + divPS +n-[P]_ = U 



(All) 



(A12) 
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where and g A are the momentum densities, P ± and 
P s , the pressure tensors and fy and f A the density force 
fields. Precise expressions for these quantities can be 
found in for the bulk and in || for the interface. 

Finally, the balance for the total energy is given by 

d 



dt 



divJ^ 



— e s A + divJ* tll 







[Je\- = 



(A13) 



(A14) 



1 d 



u s 1 u s v ^ fJ, k u 



d 



dt 



T s dt 



A s 



k=l 

dt A 



T s dt 



Pk 



A 



2( :(/ 



5« 



(A20) 



In (A19) the primes refer to the bulk baricentric velocity 
and in (A20) to the surface baricentric one. Moreover 
the thermodynamic definitions of the intensive parame- 
ters have been used. 

Following §] we will take 



where jf and J| are the energy currents. The inter- 
nal energy equations follow from the previous equations 
provided the total energy is expressed in terms of the 
internal energy and kinetic and potential contributions. 
One then arrives at 



d 

— u v + divJ± = <j± 



(A15) 



A\ = ([£>„]+, [25,,]+) 
A s 2 = ([B n ] + ,[H l{ } + ) 



(A21) 
(A22) 



as the conjugate variables of the surface polarisations. 
Then, by inserting the expressions for the balance equa- 
tions previously stated in Eqs. (A19) and (A20), one ar- 
rives at the expressions for the bulk and surface entropy 
balance equations 



dt 



divJ, 



u.\\ 



[Ju]- 



(A16) 



d 

—s v + divJf = of 



(A23) 



where J^f and J* are the internal energy current and 
± and a s u are the internal energy productions (see 0] 



o 
and 



for explicit expressions of these quantities). 



d_ 

df 



divJf 



[Js]- 



(A24) 



In order for the system of balance equations to become 
a complete set of differential equations for the fields, we 
need explicit expressions for the currents appearing in 
them. These can be obtained in the framework of non- 
equilibrium thermodynamics. 

In fact, after assuming local equilibrium, wc then pro- 
pose the fundamental equations in the bulks and at the 
surface, 



where jf and Jf are the entropy currents and of and 
of are the bulk and surface entropy productions. The 
explicit calculations for jf and J|, and of and erf may 
be found in |i| and ||. For the drift-diffusion model we 
are considering, which in particular neglects relaxation 
polarisation phenomena, generation-recombination pro- 
cesses, thermal, magnetic and viscousus effects, one then 
obtains [01 



s A = f s (u A ,plP A ,M s A ) 



(A17) 



(A18) 



(or Gibbs equations when they are expressed in differen- 
tial form) to hold locally. After derivation with respect 
to time, one then obtains the following equations in the 
bulks 



and 



1 - / 
-ji^2 J k ■ (gradfif - ZkE*) 



(A25) 



k=l 



Jf 




(A26) 



k=l 



d_ 



1 d 



and at the interface 



3 ± a 
2^ T ± dt Pk 

k=l 



dt 



P 



dt 



±' 



where the primed fluxes are defined as = 



v ± d ± 
p ir 



Here is the baricentric velocity defined through 
= Jf. For a one-dimensional system, Eq. ( |A25| ) is 
(A19) precisely Eq. tfflfy used in Section III.B. 

A similar, but much more complicated, analysis pro- 
vides the expressions for the surface flux and entropy 
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production For a one-dimensional system no equa- 
tion for the flux is necessary. For the entropy production, 
in the drift-diffusion approximation, one obtains ||, 

2 2 



k=l 



k=l 



where terms proportional to the third power of velocity 
have been neglected, and where is the surface variable 
corresponding to the normal component of the electric 
field. In this paper we have assumed that the surface 
is not polarizable, that is PJJ = 0. Due to the rela- 
tion, see p|, P% n = — E„, one then has E^ = 0. In this 



case, Eq. (A27) transforms into Eq. (J33|) used in Section 



III.B. In some cases it is useful to take into account the 
surface polarisation, po|| . In these cases, one can show 
that the discontinuity in the electric potential is given 
by [V]~ = —E^, and that one needs an additional sur- 
face state equation, that can be written, for instance as 
P%,n — P A,n(\p ]+)- K is worth noting that even in this [13] 



case Eq. (|A27[ ) can be rewritten as Eq. (pq). 
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